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Abstract 

The density profile of a hot nuclear system produced in intermediate energy heavy ion collisions 
is studied in a microcanonical formulation with a momentum and density dependent finite range 
interaction. The caloric curve and the density evolution with excitation are calculated for a number 
of systems for the equilibrium mononuclear configuration; they compare favorably with the recent 
experimental data. The studied density fluctuations are seen to build up rapidly beyond an exci- 
tation energy of ~ 8 MeV/u indicating the instability of the system towards nuclear disassembly. 
Explicit introduction of deformation in the expansion path of the heated nucleus, however, shows 
that the system might fragment even earlier. We also explore the effects of the nuclear equation 
of state and of the mass and isospin asymmetry on the nuclear equilibrium configuration and the 
relevant experimental observables. 
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I. INTRODUCTION 



In intermediate energy nuclear collisions, a hot nuclear system is initially produced in 
unstable equilibrium, generally in a compressed state. The unbalanced thermal and com- 
pressional pressure push the system towards expansion, till an equilibrium configuration is 
reached. If the system is sufficiently excited (above a few MeV per nucleon), the mononu- 
clear configuration is ultimately driven towards a subnuclear density (break-up density pf) 
below which it no longer remains in the said configuration and undergoes disassembly into 
a number of fragments. 

Study of the properties of hot nuclear systems at subnormal densities is of utmost phys- 
ical interest in understanding explosive nucleosynthesis in supernova explosions [H, 0], it 
plays a decisive factor in the element composition in the heavenly bodies. Nuclear multi- 
fragmentation studies in the laboratory serve as a window to have a closer look at these 
processes; the properties of the diluted nuclear material prior to fragmentation have a key 
role in such an understanding. It further helps in mapping the nuclear equation of state 
(EoS) and also extracting, from the complex experimental data, the density dependence 
of the nuclear symmetry energy [l, 0] that are crucial inputs in detailing the supernova 
explosion and discerning the properties of neutron star matter. 

From experimental and theoretical studies, it is seen that the break-up density pj de- 
creases with increasing excitation energy. Since this density has a profile, use of a single 
value of it leaves some room for ambiguity. In experimental studies, it is understood as an 
average density. There are several experimental estimates of the break-up density. This has 
been determined from the studies of correlation functions of emitted light particles from the 
source (5|, but the results seem to depend on the choice of the pair of correlated light frag- 
ments. It raises the question on the representation of the evolutionary stage of the reaction 
process as given by the correlation data and thus there remains a latitude for considerable 
uncertainties. The systematics of the Coulomb barriers required to fit the intermediate mass 
ejectile spectra also give an estimate of the break-up densities j6|. However, in this extrac- 
tion process, it is implicitly assumed that the barrier corresponds to a binary system Q, 
contrary to the fact that the hot system disassembles into many fragments. This introduces 
doubts in the proper estimation of the break-up densities. Recently, the break-up densities 
have also been derived from the analysis of apparent level density parameters required to fit 
the measured caloric curves [8| ; the extraction, however, hinges on a suitable mix of theory 
j9| with experiment. The deduced break-up densities from different sets of experimental 
data have sometimes wide variations among them, a sound knowledge of these densities is 
thus still incomplete. 

The disassembly of a mononuclear configuration at excitations above ~ 3 — 4 MeV/u 
has been theoretically studied in a statistical framework for more than the last two decades. 
There have been a variety of these models [l(J 11], the basic underlying premise is that 



the excited nuclear system expands to a freeze-out volume and then the system undergoes 
a one-step prompt multifragmentation, the fragmentation pattern being determined by the 
available phase space. It is assumed that both thermal and chemical equilibrium are estab- 
lished among the fragments at the freeze-out density pj so that the fragmentation pattern is 
frozen out there. In most of these models, the freeze-out density is independent of the exci- 
tation energy. Different models, however, use widely varying freeze-out densities. Whereas 
in the canonical or microcanonical models [l(J ll|, the density p/ is ~ 0.12 — 0.2 po (where 
p is the nuclear saturation density), the corresponding quantity in the lattice-gas model 
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12, ll3| is ~ 0.3 — 0.4po- This wide uncertainty in both the theoretical and experimental 



arena has prompted in recent times to have a closer look on the dependence of the nuclear 



density with excitation [14], [15|. These calculations depend on the experimental premise 
that the hot systems formed in nuclear collisions are isolated systems having a fixed exci- 
tation energy (microcanonical ensemble); for equilibrium, they are driven to an expanded 
configuration with the maximum entropy 

The calculations of Sobotka et al fl4] on nuclear expansion with excitation are schematic 
in nature; they employ a parametric form of the base density and then allow a self-similar 
expansion of the density to look for the maximum entropy configuration, the entropy being 



calculated in the Fermi-gas model. In the calculations reported in Ref. [15( the same physical 
picture was followed, however, a realistic effective Hamiltonian was used to calculate the base 
density profile in a Thomas-Fermi framework. This renders the base density self-consistent 
and allows one to calculate the entropy microscopically; furthermore, it offers the possibility 
of exploring the effect of the nuclear EoS by using a variety of effective interactions which, 
however, was not reported. 

In Ref. [HI, calculations were performed with the Skyrme interaction (SkM*) for only one 
system, namely, 150 Sm. For an expanding system pursuing the maximum entropy, the surface 
diffuseness is likely to play a key role [16|, a zero-range interaction like the Skyrme force may 
possibly not be suitable for studying such a density profile. It has been reported [l7| that a 
constrained expanded system in a Thomas-Fermi approach leads to numerical instabilities 
and the gradient (surface) terms in the energy density functional were replaced with a 
suitable Yukawa interaction [l8[ . In the present communication, we report calculations using 
the modified Seyler-Blanchard (SBM) interaction [l9, 2(|; it is of finite range and momentum 
and density dependent. Different sets of SBM parameters yield different nuclear EoS. In the 
present communication, we have explored the role of nuclear EoS on the observables related 
to an expanded mononuclear equilibrium configuration. We have further studied the isospin 
and mass dependence of the observables. In addition to expansion, a possible deformation 
along the expansion path might contribute additional entropy and mimics a fragmentation 
channel. We have considered volume-conserving quadrupole deformation to explore this 
aspect. 



II. MODEL OVERVIEW 



The key tenets of the model have been reported in brief in Ref. [15f . In the following, we 
present a somewhat more comprehensive account of the theoretical framework used in the 
present calculation. 



A. Basic Methodology 

In the experimental conditions, when two nuclei collide at intermediate energy, a hot 
nuclear system with a fixed number of neutrons and protons is formed in a nonequilibrium 
state. The system might also be compressed initially resulting in a collective flow in the 
decompression stage. We ignore the collective flow in the present calculation. As the system 
remains in isolation, the total excitation energy does not change in the subsequent evolution 
of the system, i.e., the system is microcanonical. To attain equilibrium, this hot system 
evolves in quest of maximum entropy. It is, however, still possible to describe the system 



3 



statistically by an effective temperature T. To keep the vestiges of a canonical temperature 
T in the calculation has the operational advantage that it helps in defining an occupation 
function that can be employed in evaluating various observables like energy, entropy, etc. 
To obtain the equilibrium state, we have adopted the following calculational procedure for 
a given excitation energy E*\ 

i) The system is prepared in a finite temperature Thomas-Fermi (FTTF) method such 
that the thermal excitation energy E t ^ er (T) = E(T) — E(T = 0) is less than E*. Here E(T) 
is the total energy of this system at the temperature T. 

ii) Keeping the temperature fixed at T, the system is allowed to undergo expansion till 
the total excitation (thermal-plus-expansion) of the system equals E*. The entropy for this 
configuration is evaluated. 

iii) We repeat the above steps with different initial temperatures T and from the plethora 
of different configurations so obtained, the one with the maximum entropy represents the 
equilibrium mononuclear configuration. 

In essence, this procedure aims at maximizing the entropy with respect to the collective 
coordinate that describes the expansion of the hot system under the constraint of constant 
excitation energy. Along the process there is a transfer from the initial thermal energy of 
the system to expansion energy in search of maximal entropy. 



B. Generation of the base density profile 



19, 20 



v eS (r 1 ,r 2 ,p,p) = -Ci >u 



1 - 



(1) 



The FTTF framework is employed to generate the base density profile at the temperature 
T. The effective interaction used for the calculation is taken to be the SBM interaction 
which is given by 

P 2 J2r r \, ( exp(-r/a) 
— -d {p n +p r 2 } — — — — . 
cr J {r/a) 

Here r = |rj — r 2 1 and p = |pi — P2I are the relative separation of the interacting nucleons in 
configuration and momentum space, p^x) and p(r 2 ) are the densities at the sites of these two 
nucleons, and C\ and C u are the strengths for like pair (n — n, p — p) and unlike pair (n — p) 
interactions, respectively. Following Refs. ljj 2(j, the six parameters Ci, C u , a, b, d and n are 
determined by reproducing the volume energy per particle of symmetric nuclear matter, its 
saturation density, volume asymmetry energy, surface energy, the energy dependence of the 
real part of the nucleon-nucleus optical potential and energies of Isoscalar Giant Monopole 
Resonances (ISGMR). The density exponent n governs the nuclear EoS. To explore the effect 
of EoS on the relevant observables, we have chosen two sets of SBM parameters (given in 
Table I) with n =1/6 and 4/3 that yield values of the nuclear incompressibility K OQ = 238 
and 380 MeV in infinite matter, respectively. The above interaction with n =1/6 reproduces 
quite well the ground-state binding energies, root mean square charge radii, and the ISGMR 
energies for a host of even-even nuclei varying from 16 O to very heavy systems. It has 
also been seen that for symmetric nuclear matter, the results obtained with this interaction 
agree very well (2ll . 22| with those calculated microscopically with a realistic interaction in 
a variational approach 0, [H| . 

The occupation function n T (r, p, T) for the finite system is obtained [20j from the mini- 
mization of the thermodynamic potential 

(2) 



tt = E -TS -J2^tN t . 
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It is given by 



n T (r,p,T) 



1 + exp 



2m 



+ V?(r) + p 2 V T \r) + V T \r) + 6 T , z V c (r) - pt T /T 



(3) 



In the above, E and S are the total energy and entropy of the nucleus, N T refers to the 
number of neutrons or protons (r is the isospin index), \i T their chemical potentials, and m is 
their mass, taken to be the same for neutrons and protons. The quantities Vf? and are the 
components of the single-particle potentials corresponding to the momentum independent 
and momentum dependent parts of the interaction; V 2 is the rearrangement potential arising 
out of the density dependent part of the interaction and V c is the Coulomb potential which 
contains both the direct and the exchange parts. Expressions for the different components 
of the single-particle potential are given in Appendix-A. The momentum dependent part V* 
determines the nucleon A;— mass m Tj k(r) as 



V 



+ p 2 V T \r). 



2m Tjfc (r) 2m 

The effective single-particle (SP) potential V T (r) is given by 

V T (r) = V T (r) + V T \r) + 5 T>z V c (r). 
The occupation function given by Eq. (3) is then rewritten as 

^2 



n r (r,p,T) 



1 + exp 



2m Tifc (r) 



+ V T (r) -n T )/T 



The base density is obtained from 

p T {r) = — J n T (r, p, T)dp = A* T (r) J 1/2 (rj T (r)) , 



where 



AW) = % [2m T , k (r)Tf 2 



and Jk(Vt) is the Fermi integral 

POO 



X 



K 



+ exp(x - rj T ) 



dx, 



with the fugacity i] T given as 



nr(r) = \fir-V T (r)]/T. 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



The occupation probability n T (r, p, T) given by Eq. (6) can also be expressed as the SP 
occupancy in the energy space as 



f(e T , f i T ,T) = [l + exp{(e T - f i T )/T}] \ 



(11) 
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Pr{r) 



2tt 2 



3 

/ y/e T -V T (r)f(£ T ,fj, T ,T)d£ T . (13) 

JVM) 



with the SP energy written as 

2 

e T = - P . r + V T (r). (12) 
2m Tjfc (r) 

The base density from Eq. (7) can be recast, in terms of the SP potential and energy, as 

'2m r , fc (r_ 

'V T (r) 

This density profile is generated self-consistently in an iterative procedure. The chemical 
potentials \x r are determined from the nucleon number conservation: 

N T = J p T (r)dr 

= J g T {e r ,T)f{e T ,ft T ,T)de r , (14) 
where the single-particle level density g T (e T ,T) is given by 

9r(e T ,T) = ^§J (m T>k (r))%y/e T -V T (r)r 2 dr. (15) 

The modeling of a hot finite nuclear system poses some problems as it is thermodynam- 
ically unstable. The main difficulty arises in taking proper account of the continuum states 
which are occupied at a finite temperature as a result of which the particle density does 
not vanish at large distances. The extracted observables then depend on the size of the 
box in which the calculations are performed. This problem is overcome in the so-called 



subtraction procedure [25|, [26j where the hot nucleus, assumed to be a thermalized system 
in equilibrium with a surrounding gas representing evaporated nucleons, is separated from 
its embedding environment. The method is based on the existence of two solutions to the 
finite temperature Thomas-Fermi equations, one corresponding to the liquid phase with the 
surrounding gas (Ig) and the other corresponding to the gas (g) phase. These two solutions 
are obtained from the variational equations 



and 

5Q 
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(16) 



0, (17) 



where Qi g and fl g are the thermodynamic potentials [as defined in Eq. (2)] of the respective 
systems with the same chemical potentials. The base density profile for the hot nucleus in 
question is given by p T = p T j g — p Tj9 (which may also be called the liquid profile); this is 
independent of the box size in which calculations are done. The density p T goes to zero at 
large distances implying vanishing surface pressure. 

The expressions for pi g (r) and p g (r) when written in terms of the SP potential and 
energy are analogous to that given in Eq. (13). Obviously, the nucleon k— mass and the SP 
potentials in the two phases are different and 

[Pr,i g {r) - p T , g {r)\ dr = N T . (18) 
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C. Expansion of the system 

The density profile of an expanded system, in principle, can be treated in a constrained 
Thomas- Fermi procedure [l7j]. However, for large expansions, instability sets in which is 
amplified in the subtraction procedure for the extraction of the bloated nuclear density. We, 
therefore, simulate the expansion of the system through a self-similar scaling approximation 
for the density: 

p x (r) = A 3 p(Ar), (19) 

where the scaling parameter A is unity for the unbloated nucleus and decreases with expan- 
sion, lying in the range < A < 1; p\{r) is the scaled density and p(r) is the base density 
profile generated in the FTTF framework. 

Besides its simplicity, there is no a priori reason for choosing the self-similar expansion. 
However, some justification can be found on the grounds that in a simplistic situation with a 
harmonic oscillator potential at small temperatures, the scaled and the constrained density 
profiles are equivalent. This is shown in Appendix-B. 



D. Effect of correlations on the density 

The self-consistent Thomas-Fermi procedure generates the single-particle (mean field) 
potential; it does not include the coupling of the single-particle motion with the collective 
degrees of freedom [271 ] . This coupling introduces an extra energy dependence in the effective 
mass of the nucleons (in addition to the contribution from the momentum dependence in 
the effective interaction). Incorporating this effect, the total effective mass m* (the effective 
mass is position-dependent; now on, this dependence is not explicitly shown) is then defined 



as 



m 



m ■ 



m m 



(20) 



where m w (the uo— mass) is the energy or frequency depende nt p art of the effective mass. The 
uj— mass is surface peaked and has values generally larger [28( than the nucleon mass. This 
increased effective mass has the effect of bringing down the the excited states from higher 
energy to lower energy near the Fermi surface, thus increasing the many-body density of 
states at low excitations. The system can then accommodate comparatively more entropy 
at a given excitation energy; in the present case, in searching for the state with maximum 
entropy, these correlations might have a significant role to play. An ab initio determination 
of the u— mass [28|, 2^] as a function of temperature is very involved, we have therefore taken 
a phenomeno logical form 30, 31, 32| for m u : 



^ = 1-OAAiexp 

m 



T 



21A- 



1 dp{r) 
p(0) dr 



(21) 



where p(0) is the central density of the density distribution, the temperature T is measured 
in MeV, and m is the nucleon mass. The collectivity as introduced refers to the liquid phase 
only; the density in the above equation is then p(r) = p lg (r) — p 9 (r) and A refers to the 
liquid mass. We have thus taken the a;-mass of the liquid-plus-gas phase as m 1 ^ = m u , the 
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cu-mass of the liquid; for the gas phase, the frequency dependence in the effective mass is 
neglected, i.e., m% = m. 

A self-consistent calculation of the density profile with the inclusion of u— mass is ex- 
tremely complex; to avoid the complexities, we have, therefore, adopted the following ap- 
proach which is a realistic extension of the method given in Ref. j3lj. For each fixed 
temperature T, we assume that the nucleus is well described by the mean field and the effec- 
tive mass given by Eq. (20). The self-energy of a particle in the nuclear medium is usually 
approximated by a local field (m* jm) U(r), where U(r) is a local potential (for instance, 
a Woods-Saxon potential). In our case, m^ is explicitly included in the self consistently 
generated V(r). Therefore, due to the frequency dependence in the effective mass, g(e) as 
defined in Eq. (15), is modified, taking into account the subtraction procedure, as 



AV2 

7i h 3 



m 



i g / i g( . m 

m r,k \ £r-Vr{r) 



r dr. 



(22) 



In Eq. (22), the first term in the square bracket corresponds to the liquid-plus-gas part and 
the second term corresponds to the gas part. 

The chemical potential fi T is now modified to jl T to conserve the particle number: 



N T = J g T (s T ,T)f(£ T ,fl T ,T)de T . 



(23) 



The expressions for the densities are also accordingly modified: 

3/2 



2n 2 fr? 



m 



m 

s T - V l T —f{e T , fi T , T)de T , 
ml 



(24) 



where i refers to Ig or g. 

In the present context, the following prescription is adopted to include the effect of 
correlations on the scaled density p\; the single-particle potential Va and the k— mass m fc A are 
evaluated for the scaled density, the u— mass m w A is calculated from this density according 
to Eq. (21), and then the single-particle level density g\ and the chemical potential jl\ are 
generated by using Eqs. (22) and (23) as already stated. This gives the correlated scaled 
density p\ T {r). 



E. Energy and entropy 

In absence of correlations, the total energy of the finite unexpanded nucleus is expressed 
in terms of the uncorrelated density and the occupancy given by Eqs. (13) and (6). The 
total energy density is obtained as 

^l^PW + ^W], (25) 
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where E^ in (r) and £J. nt (r) are the kinetic energy and the interaction energy density, respec- 
tively The kinetic energy density is 

^2 



2tt 
mh 3 



(2m T , k Tf 2 J 3/2 Mr)), 



(26) 



where the effective mass m Tjfc and the Fermi integral Jr(v) are defined through Eqs. (4) and 
(9), respectively. 

The interaction energy density is given by 



V T °(r)p T (r) + ^V T \r) J p 2 n T (r,p)di 



+ £c(r)5 TiZ 



= \v?{r)Pr{r) + ^ (2m T , k Tf 2 J 3/2 (r, T (r)) V T \r) + E c (r)6 T , z . (27) 



The Coulomb energy density E c (r) is the sum of the direct and exchange contributions: 



The direct term is obtained as 



£f(r) = Tre 2 ^^ f°° Pp (r') [(r + r') - |r - r'|] r'dr 7 , 
r Jo 



and the exchange term is calculated from the Slater approximation: 



^p 2 

__( 3vr 2^1/3 n 4/3 
47T 



Combining Eqs. (26) and (27), one may write (25) as 



2ir 



(2m T , k Tf 2 J 3/2 ( V r(r)) i— + V*)+ ^V T °(r)p T (r) + E c (r)& 



m 



The total energy is 



E 



£(r)dr. 



(28) 
(29) 

(30) 
(31) 

(32) 



The total entropy in the Landau quasi-particle approximation for the nuclear system with 
the base density profile is 



S = -Hf 9r(e T ,T) [/In /+(!-/) In(l - /)] de T , 



(33) 



where the SP level density and occupancy are given by Eqs. (15) and (11). 

The energy of the expanded correlated system is calculated with the suitably modified 
effective mass, single-particle potentials and fugacity. The entropy of the correlated system 
can be similarly calculated. In the subtraction procedure, the energy E of the nuclear liquid 
is obtained through 

E = E lg - E g , (34) 

where E ig and E g are the total energies of the liquid-plus-gas system and of the gas alone; 
these energies are calculated from Eq. (31) with appropriate modification for expansion and 
correlation effects. For the entropy, the subtracted SP level density of the form given by 
Eq. (22) (that includes correlation effects), but modified for scaling is employed. 
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F. Deformation of the system 

So far, we have taken the shape of the hot expanding system to be spherical. In search of 
the maximum entropy configuration, the system may, however, explore possible deformation 
paths along with expansion mimicking fragmentation. For simplicity, we consider only a 
volume conserving quadrupole (spheroidal) deformation at all stages of expansion and check 
whether there is additional gain in entropy from deformation of the system. 

In a volume conserving deformation, only the surface and Coulomb energies change. To 
calculate these changes, a sharp surface approximation to the density profile is made with 

-Rsharp — y |(r 2 ). This also facilitates the calculation of the entropy gain from deformation. 

The excess Coulomb energy 5E c (X,j3) of the expanded deformed system is given by [3" 



5E C (\,(3)=E C (\,0) [(1 



x. 



1/31^1-1/2 



F{x) - 1] , 



(35) 



where /3, the deformation of the spheroidal system, in terms of the lengths 2a and 2c of its 
principal axes, is defined as 



0=3 ln(c/a), 



(36) 



where a < c for prolate and a > c for oblate spheroids. The constraint of volume conservation 
gives a 2 c = -Rs harp . The quantity x, in terms of /3, takes the form 



3 (1 + /3/8) 
2 1 (1 + /3/2) 2 



(37) 



and the function F(x) reads 



1 



F(x) = - [ln(l + y/x) - ln(l - y/x)] , for a < c, 



tan 



x 



1/2 



for a > c. 



(38) 



The expression of the surface area of the deformed nucleus is 

2irac 2 



A(X,(3) 



2na z 



sin 



for a < c, 



7T 



y/a 2 — c 2 



2a 2 y/ a 2 — c 2 + ac 2 In 



+ yfc 



,a — y/a 2 — c? 

The change in surface free energy due to deformation is given by 

6F(\,/3) = 6A(\,I3)<t(p,X,T), 



for a> c. (39) 



(40) 



where 5 A is the excess surface area due to deformation. The surface tension coefficient a is 
density, temperature, and asymmetry dependent; it is taken as 



a(p, X, T) = a{p , X,T — 0) y{p) a{T), 



(41) 



with a(po,X, T = 0) = cr(p ,X = 0,T = 0) — a s X 2 ; here a s is the surface asymmetry 
coefficient taken to be 1.7826 MeVfm -2 , X = (N — Z)/A is the asymmetry parameter of 
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the nucleus, and a(po,X = 0, T = 0) = 1.06 MeVfm -2 is the surface energy coefficient of 
symmetric semi-infinite nuclear matter. 

The surface tension has its maximum value at the saturation density po- The density 
dependence of the surface tension coefficient is evaluated by applying the scaling approxi- 
mation on the ground-state density profile of semi-infinite nuclear matter. In the sub-nuclear 
density region, it is found to be well approximated by a polynomial y(p) of the form 

4 

y{p) = Y. h ^p/p») k - ( 42 ) 

k=l 

The values of the b k coefficients are h x = 0.5382, b 2 = 2.3124, b 3 = -2.2312, and b 4 = 0.3816. 
The temperature dependence of the surface tension is taken as 35 

^HIK'-If ■ (43) 

with T c = 15 MeV, the critical temperature for nuclear matter. 

The excess entropy from deformation is calculated as 5S = —d(5F)/dT\ p . The excess 
surface energy is obtained from 5E snT f = 5F + T5S. The total energy of the expanded 
deformed system can then be written as 

E(X, (3, T) = E(X, 0, T) + SE C (X, (3, T) + 5E sur{ (\, (3, T). (44) 

The excitation energy is given by 

E*(X, (3, T) = E(X, (3, T) - E{X = 1, (3 = 0, T = 0) (45) 

and the corresponding total entropy is 

S(X, (3, T) = S(X, 0, T) + 5S(X, (3, T). (46) 
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III. RESULTS AND DISCUSSIONS 



To explore the density reorganization in nuclei with increasing excitation energy, we have 
chosen three isobars of mass v4=150 (namely, Cs, Sm, and Yb) and two isobars of mass 
^4=40 (Ca and S) as representative systems. The two different sets of isobars with different 
isospins are chosen in order to study the mass as well as the asymmetry dependence on the 
evolution of density and some other characteristic observables related to the mononuclear 
configuration. Two sets of SBM interaction parameters (given in Table I) with widely varying 
nuclear incompressibility are used to investigate the effect of nuclear EoS on the observables 
considered. All the calculations reported are performed at different fixed excitation energies 
with inclusion of thermal and expansion effects and also with incorporation of the correlation 
effects introduced through the frequency dependence of the nucleon effective mass. 

The dependence of the observables related to the excited mononuclear configuration on 
the nuclear EoS has not been studied so far. In Fig. 1, the density evolution of the hot 
Sm nucleus is displayed as a function of the excitation energy E*/A for two different EoS, 
one softer with = 238 MeV and the other one harder with = 380 MeV. Here, 
the quantity p c represents the central density of the hot expanded nucleus in equilibrium 
corresponding to the maximum entropy configuration and po is the ground-state central den- 
sity evaluated in the self-consistent Thomas-Fermi framework. Obviously, at zero excitation 
Pel Po — 1- It is found that for the softer EoS, above E*/A ~ 11 MeV, the entropy increases 
monotonically with expansion and hence the mononuclear configuration ceases to exist; for 
the harder EoS this limiting excitation energy is somewhat higher. At a given excitation, 
the equilibrium density is seen to be lower for the softer EoS. This can be understood from 
the following qualitative argument: the maximum entropy configuration is an outcome of 
the balance between the entropy gain from expansion and a loss due to reduction in temper- 
ature. Since for a harder EoS, the gain in entropy from expansion is energetically costlier, 
for a given excitation energy, the expansion is less for a harder EoS. A more quantitative 
explanation in a simplistic model is given in Appendix-C. The effects of correlation on the 
mononuclear observables with increasing excitation energy for a Skyrme-type (SkM*) force 
have been discussed in an earlier paper [15]. It has been found that the changes in the 



relevant observables with the inclusion of correlation effects with the SBM force are quali- 
tatively the same as the SkM* interaction and therefore we do not want to delve into these 
details. However, for completeness, we present some selective results only for the nucleus 
150 Sm. In Fig. 1, the dashed curve displays the results for the evolution of central density 



without correlation effects for the softer SBM force. As in Ref . [15( , the effect of correlation 
is found to reduce the central density at all excitations considered. 

Two sets of experimentally derived data are also shown in Fig. 1; the filled circles are the 
experimental points extracted from the apparent level density parameters [8] for the mass 
selection 140 < A < 180 and the open squares are the ones obtained from the Coulomb 
barrier systematics [6j for Au-like systems. It is found that the calculated results with the 
softer EoS having nuclear incompressibility very close to the generally accepted value 36| 



compare favorably well with those derived from the level density parameters. 

In the following, the calculations are done with the softer EoS and the reported results 
are to be understood in its terms, unless otherwise mentioned. In the upper panel of Fig. 2, 
the density (p c ) evolution with excitation for the three nuclei Cs, Sm, and Yb is shown 
along with the experimental data 0,0]. Since the density has a profile, its characterization 
by a single entity p c leaves some room for ambiguity. We have, therefore, also considered 
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an average density p = A/V eq where V eq is the effective volume of the excited nucleus in 
equilibrium, taken as V eq = ^B? sharp where the sharp surface radius R s har P is calculated as 

Rsharp = y |(r 2 ) from the density distribution. In the lower panel of the figure, the density 

ratios are displayed as Vo/V eg as a function of excitation energy per nucleon. Here Vq is 
the ground-state effective volume. Inspection of the two panels shows that the two sets of 
density ratios so defined are somewhat different. The difference is more manifest at higher 
excitations. The experimental data obtained from the apparent level density parameters lie 
in between the two sets of calculated densities. 

Fig. 2 displays the asymmetry dependence of the density evolution with excitation. It 
is seen that the more symmetric nucleus has generally a higher equilibrium density at a 
given excitation. This may be understood by examining the asymmetry dependence of the 
incompressibility Ka of finite nuclei. It can be written in powers of A^ 1 ^ 3 with a liquid-drop 
type expansion (36L I37L l38"|: 

K A = K V + K S A~ 1 ' 3 + K 5 5 2 + K C Z 2 A~^ 3 + ... . (47) 

Here K v , K s , Kg and K c are the volume, surface, asymmetry and Coulomb coefficients, 
respectively, and S is the asymmetry parameter 5 = (N — Z)/A. For isobars, the asymmetry 
and Coulomb terms are different. From the known [37], [38[ values of Kg and K c (which are 
negative) , it is found that the more symmetric nucleus has a larger incompressibility Ka and 
thus, in the context of the discussions in relation to Fig. 1, it follows that the equilibrium 
density of the more symmetric nucleus is larger. For the nucleus with the largest charge 
considered here, there is a marked deviation from the smooth decrease of density with 
excitation in the energy range 6 < E*/A < 8 MeV; it is difficult to discern the exact reason 
for it because of the delicate interplay of the Coulomb, thermal and expansion effects on 
the density distribution. The manifestly opposite behavior of the density ratios as defined 
in the two panels of Fig. 2 can be understood by examining the density distributions at the 
relevant excitations as given in Fig 3. Particularly noticeable is the inversion of the trend 
of the central density with excitation at E*/A = 7 MeV. The steep decrease of the density 
ratio given by Vo/V eq can be understood from the rapid growth of the density tail with 
excitation. The evolution of density with excitation energy for the A = 40 systems (Ca and 
S) is displayed in Fig. 4. As in the case of heavier systems explored (Fig. 2), the density 
for the more symmetric system is higher; the asymmetry effect on the density is seen to be 
more prominent for the heavier systems. It is further found that for the same asymmetry, 
the heavier systems have somewhat higher equilibrium density at a given excitation. It is 
noted that the two definitions used for the characterization of density yield values which are 
not much different for the lighter systems. 

The expansion of the nucleus has an important bearing on the correlation of excitation 
energy with temperature. The correlation (caloric curve) so obtained for the expanded hot 
nucleus 150 Sm is displayed in Fig. 5. At a fixed excitation energy, the system cools down 
with expansion, and therefore the recorded temperature at the equilibrium configuration is 
significantly lower than that for the unexpanded nucleus prepared initially with the same 
excitation. In all of our calculations, the temperature refers to the canonical temperature. 
We have checked that the microcanonical temperature obtained from T~ 1 = dS eq /dE*, is 
not much different from the canonical one. To explore the role of the EoS on the correlation 
between E* and T, the caloric curves obtained from both the softer (K^ = 238 MeV) and 
the harder (K^ = 380 MeV) EoS are shown in Fig. 5. It is seen that at a given E*, the 
equilibrium temperature is significantly higher, particularly at high excitations, for the stiffer 
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EoS. Tentatively, this can be understood from the fact that the equilibrium density is higher 
for a harder EoS (Fig. 1), the corresponding expansion is less and so more energy is locked 
in the thermal mode with a resultant higher temperature. The same conclusion is arrived at 
from the simplistic Fermi-gas type framework presented in Appendix-C. A representative set 



of experimental data [39| for medium-mass nuclei are also shown in the figure. It is seen that 
the calculated caloric curve with the softer EoS compares very well with the experimental 
data. The dashed curve represents the caloric curve without inclusion of correlation effects 



{m^/m = 1). As in Ref. [15J], introduction of correlation is seen to reduce the temperature 
for the equilibrium configuration at all excitation energies considered. 

The dependence of the nuclear caloric curve on mass and isospin asymmetry is shown in 
Fig. 6. The top panel refers to the heavier isobars and the bottom panel to the lighter ones. 
It is apparent from the figures that the caloric curves are nearly independent of asymmetry 
at lower excitations. For E*/A > 5 MeV, the isospin dependence becomes increasingly 
prominent with increasing excitation, the equilibrium temperature is higher for the more 
symmetric system. The mass dependence on the caloric curve is seen from an examination 
of the results for pairs of nuclei which have almost the same asymmetry like Sm and S or Yb 
and Ca. It is found that the equilibrium temperature for the lighter nuclei are comparatively 
always lower at all excitations. This behavior of the caloric curve with mass and isospin can 
be explained from the fact [using Eq. (47)] that for the same asymmetry, heavier nuclei have 
larger incompressibility; also, for isobars, the ones with more asymmetry have an effectively 
lower incompressibility. 

The caloric curves show an interesting feature, namely, the occurrence of negative heat 
capacity, generally beyond an excitation energy E*/A ~ 8 MeV. Intuitively, one understands 
that if a system with a given excitation expands, it does so at the cost of thermal energy 
and hence there may be a density region when temperature may decrease with increasing 
excitation if the system expands much in pursuit of maximum entropy. In a simplistic model 
in Appendix-C, we show that this occurs in the density region 1/4 < p/po < 5/8. It is also 
found that the more symmetric systems have lesser bends in the caloric curves (for Yb, in 
the excitation energy range explored, the bend shows up only later at E*/A ~11 MeV), 
this is because they offer more resistance towards volume expansion as already explained. 

The expansion energy -E exp n constitutes a part of the total excitation. It is defined as 
-Eexpn = E(X, T) — E(X = 1, T). In Fig. 7, it is displayed as a function of E* for the different 
systems. In the bottom panel, the effect of the EoS on the expansion energy is shown for 
the nucleus 150 Sm. It is seen that E expn decreases with stiffness of the EoS; the same is 
true for all other systems studied. The filled circle refers to an experimental estimate 0] 
for medium-mass nuclei (A ~ 160 — 180). This is very close to our calculations with the 
softer EoS for the 150 Sm nucleus. The top and middle panels refer to calculated results for 
isobars with A=150 and 40, respectively. As already stated, the more symmetric nuclei have 
effectively a harder EoS, this is reflected in the comparatively lower expansion energy for 
these nuclei. It is seen that in general, the expansion energy comprises a significant part of 
the total excitation. 

The state with maximum entropy S(X eq ,T eq ) corresponds to the state with maximum 
probability. Other configurations with different values of the scale parameter A and temper- 
ature T are also probable. Their probability is determined by the entropy profile S(X, T) at 
a fixed excitation energy: 

W(X,T) oc e 5(A ' T) . (48) 
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The nth moment of the central density p c is then 



Similarly, for the temperature T one has 

(T n (\,T))= J $eS{x : T)dT } • (50) 

Equations (49) and (50) allow one to calculate the mean and the variance of p c and T at a 
fixed excitation. Figs. 8 and 9 display the entropy profile for the systems 150 Cs and 150 Yb at 
different fixed excitations E* j A = 4, 6, 8, and 10 MeV as a function of the central density 
(Fig. 8) and temperature (Fig. 9) of the expanded configurations. The arrows in the figures 
indicate the configurations with the maximum probabilities, the full ones refer to Yb and 
the broken ones correspond to Cs. The entropy profiles are seen to flatten with excitation. 
This implies larger fluctuations at higher excitations. 

For a thermodynamic system, the average and the most probable (equilibrium) value of 
an observable are the same. For a finite system, however, they may differ. Experimentally, 
it is the average value that one measures. The average temperatures for the A=150 isobaric 
systems are displayed in panel (a) of Fig. 10; the corresponding quantities for 40 Ca and 40 S 
are shown in panel (c). Panels (b) and (d) display the variances in temperature for the 
isobaric systems with A=150 and 40, respectively. The equilibrium temperatures as shown 
in Fig. 6 are somewhat larger than the corresponding average values. Fluctuations build up 
with excitation. At a given excitation energy, the fluctuations for the lighter systems are 
comparatively higher as expected. For the heavier systems, a sudden increase in fluctuation 
is observed beyond E* / A ~ 9 MeV, whereas the increase is relatively smooth for the lighter 
systems. In Fig. 11, the average values of the specific volume v c (=l/p c ) in units of t>o 
(= 1/po) an d their variances for the systems considered are shown. The equilibrium volume 
is found to be somewhat less than the average value. As in the case of temperature, the 
fluctuations in volume rise smoothly up to an excitation energy E* j A ~ 9 MeV, beyond 
which the build-up is very sudden. This sudden build-up is more pronounced in comparison 
to that in temperature for both light and heavy systems. This large density fluctuation 
indicates that beyond E* / A ~ 9 MeV, the systems become unstable and break up into 
many pieces. It turns out that the negative branch of the heat capacity and the onset of 
large fluctuations occur at around the same excitation energy indicating a possible close 
correlation between them. 

Along with expansion, the excited nucleus may undergo deformation if that is profitable 
from the entropy considerations. From the interplay of Coulomb and surface energies, a 
barrier is found along the deformation path; the barrier decreases with increasing temper- 
ature and decreasing density (expansion). At a particular temperature, for the deforming 
system, in our calculations the scale parameter A is taken such that the total excitation at 
the top of the barrier matches the given excitation E* j A. This is repeated for different tem- 
peratures and the maximum entropy among these different configurations is selected. If this 
entropy exceeds that for the expanded spherical equilibrium configuration, then deformation 
is favored leading to the fragmentation channel. 

The extra entropy AS gained from deformation over that at spherical equilibrium shape 
is displayed in the upper panels of Fig. 12 for the nuclei Sm and Cs as a function of excitation 
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energy. As mentioned, we have considered a volume-conserving quadrupole deformation. For 
E*/A less than ~ 4 MeV, AS is negative in the restricted deformation space chosen; it is 
thus seen that only above this excitation domain, fragmentation resulting from deformation 
is more favorable. With deformation, the maximum entropy configuration occurs at a lower 
temperature at a given excitation; the corresponding caloric curves are shown in the bottom 
panels. In all the calculations, it is found that the oblate shapes are inhibited due to their 
lower entropy. 



IV. CONCLUDING REMARKS 

In the present work, the gross features of an expanding mononuclear configuration at 
moderate and higher excitations have been dealt with in a semi-microscopic framework. A 
finite range, momentum and density dependent realistic effective interaction has been em- 
ployed for our investigation. Renormalization of the nucleon mass coming from the coupling 
of the single-particle motion with the surface degrees of freedom is seen to play an important 
role in the correlation of entropy with excitation in a microcanonical formulation, this has 
been taken into account phenomenologically to avoid complexity. The density dependence 
in the interaction leaves room to vary the nuclear equation of state well within the ambits 
of accepted nuclear parameters. We have explored the influence of the EoS on mononu- 
clear observables and find that at a fixed excitation, a nucleus with a softer EoS has more 
expansion, the corresponding expansion energy is found to be larger and thus its effective 
temperature is comparatively lower. 

The resultant mononuclear caloric curve and equilibrium density compare very well with 
the experimental data with the chosen softer EoS with nuclear incompressibility A^ =238 
MeV, which is close to its well-accepted value. We have further studied the mass and isospin 
dependence on the mononuclear observables. For nuclei with same asymmetry but higher 
mass, it is found that the system expands less and consequently has a higher temperature; 
this is also true for nuclei with same mass but less asymmetry. All these observations 
could be explained from the fact that heavier and symmetric systems have effectively higher 
incompressibility. 



In the experimentally constructed nuclear caloric curves [39j, |41(, a plateau is observed in 
the excitation energy range ~ 3 — 8 MeV/ A. This is indicative of infinite specific heat. In 
some analyses, a negative specific heat [42j in this energy domain is also suggested. These 
are pointers to a possible liquid-gas phase coexistence. In our calculations, the plateau comes 
naturally from nuclear expansion, so also its bending (negative specific heat) at a relatively 
higher excitation. The fluctuations in temperature and volume build up with excitation. 
The rapid growth of these fluctuations at E*/A ~ 9 MeV is suggestive of the instability of 
the mononuclear configuration against prompt multifragmentation. The occurrence of the 
negative heat capacity and the sudden growth of fluctuations at around the same excitation 
energy indicate a close correlation between them. 

While expanding, the nuclear system may gain further entropy from fragmentation into 
a number of pieces, mostly from the generation of new surfaces. We have dealt with this 
in a simplistic way through a volume-conserving quadrupole deformation as a precursor 
to fragmentation. It is found that above E*/A ~ 4 MeV, the systems generally favor 
deformation. The evolution of the mononuclear configuration into a number of fragments 
may be treated through higher multipole deformation; it would be an involved task, but it 
is worth investigating. 
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Appendix-A 
The single-particle potentials 

The components of the single-particle potentials for the SBM interaction (1) entering in 



Eq. (3) can be written as [20 



V T °(n) = -2vra 2 / [1-d 2 {p( n ) + p(r 2 )} n ] u(r u r 2 ) [C,p T (r 2 ) + C u p_ T (r 2 )] r 2 2 dr 2 

u(r h r 2 ) Q (2m Tik (r 2 )T) 5/2 J 3/2 (ry T (r 2 )) 



+ 



irb 2 h 3 



+ C u (2m_ r , fe (r 2 )T) 5/2 J 3/2 (^_ T (r 2 )) 



r\dr 2 , 



(A-l) 



V?(n) 



2na A 



2 r°° 



u(r 1 ,r 2 ) [Cip T (r 2 ) + C u p^ T (r 2 )] r\dr 2 , 



(A-2) 



V 2 (n) = 2nna 2 d 2 J2{[CiPr'(n) + C uP ^( ri )] 

t' 

POO 

x/ Pr'( r 2)«( r i,r 2 ) [p(n) + pi^)}^ 1 rldr 2 }, (A-3) 

with 

p(r) = p T (r) + p_ T (r) and u(n, r 2 ) = [ e ~^~ r ^ a - e -(^+^)/ a ] /(r x r 2 ). (A-4) 

In these equations, if r refers to proton, — r refers to neutron and vice versa. From the 
structure of Eq. (A-3), it is seen that the rearrangement potential V 2 {r) is isospin indepen- 
dent. 

The Coulomb single-particle potential is the sum of the direct and exchange contributions: 

V e (r 1 ) = V?(r l )+V e E *(r 1 ), (A-5) 

where 

2ne 2 f°° 

V c D ( r i) = / [( r i + r 2) - \n - r 2 \] p p (r 2 )r 2 dr 2 , (A-6) 

r i J o 

and 

^ (ri) = _ e 2 (3/7r) l/3 p l/3 (ri) . (A . 7) 
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Appendix-B 
Equivalence of scaled and constrained density 

In this Appendix we show the equivalence between the scaled density and the constrained 
density for a harmonic oscillator (HO) potential. 

From Eqs. (7) and (8), the density at a finite temperature for a nucleus (ignoring the 
isospin) can be written in the form 

p(r) = ^(2mT) 3 /V 1/2 (^), (B-l) 

where r) = (p — V)/T, V being the single-particle potential. For small temperatures, i] is 
large. When p < V, 77 is large negative and then J\/2{v) ~ eV leading to vanishing density. 
When p > V, r\ is large positive and then Ji^iv) — f 7 ? 3 ^ 2 ? an d the density is given by 

^»^(f) 3/2 ^- v '» 3/2 < B - 2 > 

For a harmonic oscillator potential with oscillator frequency u, V — ^mu 2 r 2 , the corre- 
sponding density is 

/ l \ 3/2 

p(r) = B I fj, - -mu V J , (B-3) 
where B = {^r) 3 ^ ■ The total number of particles for this density is 

"< r ><H(£r (B - 4 » 

For the density given by Eq. (B-3), following Eq. (19), the scaled density is 



3/2 

Px(r) = B' 



/iA 2 - l -m{uo\ 2 fr 2 



(B-5) 



It is seen that the chemical potential and the oscillator frequency are now scaled as 

fix — > p\ 2 and u\ — > cu\ 2 . (B-6) 

The constrained Hamiltonian density constraining the particle number and the rms radius 
is given by 

H(r) = £ kh \r) + (V- ar 2 ) Pa (r) - p aPa (r), (B-7) 

where S km is the kinetic energy density, p a constrains the total number of particles and a 
constrains the rms radius to given values A and R a , respectively, so that 



= J p a (r)dr, (B-8) 
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and 



r 2 p a (r)dr. 



(B-9) 



In the absence of the constraint on the rms radius, p a — > p and i? a — > _R, the unconstrained 
radius. 

In analogy to Eq. (B-2), at low temperatures, the density corresponding to the Hamilto- 
nian given by Eq. (B-7), for a HO potential is written as 



p a (r) = B 



Pa — -muj 2 r 2 + ar 2 



3/2 



;b-io) 



which can be recast as 



p a (r) = B 



,2^2 



3/2 



;b-ii) 



where uj 2 a = u 2 — 2a/ m. Again, following Eq. (B-4), the total particle number from Eq. (B-8) 
is given by 



2 / p c 



3 \hw c 

From Eqs. (B-9) and (B-12), it then follows that 

, n2 ft 2 3 f3A\ 1/3 

This is also true for the unconstrained density (a =0): 



2m 2 V 2 y 



1/3 



From Eqs. (B-13) and (B-14), 



fuv a = 

^a 



Similarly, from Eqs. (B-12) and (B-15), 



(B-12) 



(B-13) 



;b-h) 



(B-15) 



R 2 

Pa = A* pa ■ 

Therefore, the constrained density given by Eq. (B-ll) is rewritten as 



(B-16) 



Pa{r) = B 



Rl 2 



2\ 2 



R2 1 I 

/x— m c<j— r 



i? 2 \ 



i? 2 / 



3/2 



(B-17) 



Identifying A as R/R a , one can see the equivalence of the scaled density given by Eq. (B-5) 
with the constrained density. 
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Appendix-C 

EoS dependence of equilibrium density and temperature 

In this Appendix, we show in a simplistic model why the density and temperature at 
equilibrium are higher for a stiffer EoS at a given excitation E*. 
The expansion energy for the decrease of density from p to p is 



E 



cxpn 



^(I-P/Po) 2 . 



(C-l) 



This is of the same form as suggested by Friedman Ka is the incompressibility of the 
nucleus. 

The thermal energy of a hot nucleus, in a Fermi-gas type model is 



-Ether — aT , 

where the density dependence of the level density parameter is given as 

n 2 



Ae F (p) 



The latter can be written in the form 



a = bp 



-2/3 



(C-2) 



(C-3) 



(C-4) 



where b is a constant. The total excitation energy is given by the sum of the expansion and 
the thermal energy, 



E* 



-Eexpn Ei 



ther • 



(C-5) 



The entropy is given by S(p, T) = 2aT = 2bp- 2 ^T . At equilibrium, S is a maximum. 
Because of the constraint that the sum of the thermal and expansion energies is equal to 
the given E*, the quantity to be maximized is 



S{p,T) = 2bp- 2/3 T-- 

v 



bp- 2 ' 3 T ,J 
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1 - p/p ) 2 - E* 



(C-6) 



where v is a Lagrange multiplier. Then dS/dT = yields v = T, and OS /dp = leads to 

K A 



bT l 



resulting in 



E 



ther 



6p 

_Ka£_ 
6 Po 



p 5/3 (l-p/po), 



(C-7) 



(C-8) 



From Eqs. (C-l) and (C-8), the total excitation is 

P o 



E* 



Ka 
18 



1 + 



Po 



(C-9) 
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Thus at a fixed E* , if Ka is higher, the quantity in the square bracket should be smaller. 
This is true if p/po increases with Ka, in the simplistic model this is satisfied only in the 
range 1/4 < p/po < 1, which, however, covers almost the whole density range of our concern. 
From Eq. (C-7), 



dT 



K A p 



2/3 



dE* 36p bT _ 



5-8^ 



dp 



PoJ dE*_ 



(C-10) 



Eq. (C-10) shows that dT/dE* is negative if the quantity in the square bracket is negative. 
From Eq. (C-9), one finds that dp/dE* is negative in the density range 1/4 < p/po < 1. In 
this schematic model, one thus finds that for 1/4 < p/p < 5/8, the temperature decreases 
with excitation. 
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TABLE I. The parameters of the effective interactions used in the present paper (in 
MeVfm units) and their bulk incompressibility modulus (in MeV). 



n 


a 


C u 


a 


b 


d 




1/6 


291.7 


910.6 


0.6199 


928.2 


0.879 


238 


4/3 


210.7 


653.2 


0.529 


615.4 


0.924 


380 
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Figure Captions 

Fig. 1 The calculated equilibrium density of 150 Sm as a function of excitation energy per 
nucleon. The thick line corresponds to the calculations with the soft EoS (^00=238 
MeV), and the thin line refers to those with the hard EoS (i^oo=380 MeV). The 
dashed line refers to calculations without correlation effects for the softer EoS. The 
filled circles are the experimental data from Ref. the open squares are the ones 
from Ref. 

Fig. 2 The calculated equilibrium densities as a function of excitation energy for A = 150 
isobars in the two definitions of density as explained in the text. The experimental 
points are the same as in Fig. 1. 

Fig. 3 The equilibrium proton density profile for the system 150 Yb at excitation energies of 
6, 7, and 8 MeV per nucleon. 

Fig. 4 Same as in Fig. 2 for the nuclei 40 Ca and 40 S. 

Fig. 5 The mononuclear caloric curve for the system 150 Sm. The thick and thin lines refer 
to soft and hard EoS, respectively. The dashed line corresponds to the caloric curve 
calculated without correlation effects for the softer EoS. The experimental data 
points (filled circles) are from Ref. [39l ] . 

Fig. 6 Mononuclear caloric curves showing isospin effects for A = 150 isobars (top panel) and 
A = 40 isobars (bottom panel). The experimental data are as in Fig. 5. Information 
on the mass dependence of the caloric curves can be gained by comparison of the 
curves of the pair of nuclei 150 Sm and 40 S that have asymmetry X ~ 0.2, with those 
of the nuclei 40 Ca and 150 Yb that are symmetric or nearly symmetric. 

Fig. 7 Expansion energy as a function of total excitation energy for A = 150 isobars (top 
panel) and A = 40 isobars (middle panel). The EoS dependence of the expansion 
energy for 150 Sm is illustrated in the bottom panel. The filled circle corresponds to an 
experimental estimate 0. 

Fig. 8 Entropy profile at different fixed excitations as a function of the central density for 
the nuclei 150 Yb and 150 Cs. The equilibrium configurations are marked by arrows, full 
ones for Yb and the broken ones for Cs. 

Fig. 9 Entropy profile at different fixed excitations as a function of temperature for the 
systems 150 Yb and 150 Cs. The arrows have the same meaning as in Fig. 8. 

Fig. 10 The average values of temperature as a function of excitation energy are shown in 
panels (a) and (c) for A = 150 and A = 40 isobars, respectively. In panels (b) and 
(d), the corresponding variances in temperature are shown. 

Fig. 11 The average values of the specific volumes and their variances for the different systems 
considered are displayed. 

Fig. 12 The entropy gain AS from deformation as a function of E*/A for the systems 150 Sm 
and 150 Cs are displayed in the upper panels. The caloric curves with (full line) and 
without (dashed line) deformation are shown in the lower panels. 
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